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Abstract 

We experimentally investigate the nonlinear response of a multilayer graphene resonator using a 
superconducting microwave cavity to detect its motion. The radiation pressure force is used to drive 
the mechanical resonator in an optomechanically induced transparency configuration. By varying 
the amplitudes of drive and probe tones, the mechanical resonator can be brought into a nonlinear 
limit. Using the calibration of the optomechanical coupling, we quantify the mechanical Duffing 
nonlinearity. By increasing the drive force, we observe a decrease in the mechanical dissipation 
rate at large amplitudes, suggesting a negative nonlinear damping mechanism in the graphene 
resonator. Increasing the optomechanical backaction, we observe a nonlinear regime not described 
by a Duffing response that includes new instabilities of the mechanical response. 
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The unique properties of graphene such as atomic thickness, low mass density, and high 
modulus of rigidity make it very attractive material for nanoscale electromechanical systems 
(NEMS) for several technological applications. After the first demonstration of few layer 
thick graphene NEM^, there has been an extensive studies on graphene nanoelectromechan¬ 
ical systems ranging from electromechanical resonatord^El^ oscillator^ and optomechanical 
systems aiming to probe the quantum regime of graphene motionP®. In this pursuit, large 
mechanical quality factors in graphene based NEMS have been demonstrated as welP®. Due 
to its atomic thickness, graphene based NEMS also exhibit rich nonlinearity such as onset of 
Duffing nonlinearity and nonlinear damping at realativly small mechanical amplitude^^*^. 
These properties further makes graphene an attractive candidate for developing optome¬ 
chanical systems to reach the quantum regime of graphene motiorpl, to store microwave 
photon^, and could possibly be useful to understand dissipation in graphene NEMS for 
improved device performance^. 

The coupling between mechanical resonator and optical/superconducting microwave cav¬ 
ities has enabled the detection of mechanical motion with excellent sensitivities^®®, offering 
an attractive platform to characterize the nonlinear response of mechanical resonators. In 
this Letter, we study non-linear dynamics of a multilayer graphene resonator by means of 
coupling it to a superconducting microwave cavity. The graphene resonator is driven by 
injecting two microwave tones in the cavity, which are detuned by the mechanical resonant 
frequency leading to an oscillating radiation pressure force which drives the mechanical res¬ 
onator. By changing the amplitude of these tones, we can independently control the driving 
force and dissipation due to the optomechanical backaction forces. We drive the mechanical 
resonator into the Duffing regime and characterize the nonlinearity. With increase in the 
driving force, we observe a reduction in linear dissipation rate, large hysteresis with sweep 
direction, and an instability in the mechanical amplitude. 

Our device consists of a multilayer graphene resonator coupled to a superconducting 
microwave cavity as studied previouslj^. Fig. la shows a scanning electron microscope 
image of a multilayer graphene resonator coupled to a superconducting microwave cavity. 
The multilayer graphene mechanical resonator is 10 nm thick and is suspended above a 
gate electrode of the microwave feedline by approximately 150 nm. The superconducting 
cavity is in a quarter wavelength coplanar waveguide geometry fabricated with an alloy of 
molybdenum and rhenium (T^ 9 K) on an intrinsic silicon substratd®. The measurements 
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are performed in a dilution refrigerator under vacuum at 14 mK. The superconducting 
cavity has a resonance frequency of cUc = 27r x 5.90054 GHz, with an internal dissipation rate 
Ki = 2n X 54 kHz and coupled to a feedline with an external coupling rate Kg = 27r x 188 kHz 
(coupling fraction r] = ^ = 0.78). The graphene resonator forms a mechanically 

compliant capacitor to the microwave feedline as shown schematically in £g. lb. Motion of 
graphene resonator modulates the capacitance and hence the cavity frequency. The graphene 
resonator has a resonance frequency of Um = 27r x 36.233 MHz. Using thermal noise, we 
calibrate the optomechanical coupling defined as go = where Xzpf are the quantum 

zero-point fluctuations of the mechanical resonator to be 27r x 0.83 Hd^, which also provides 
an absolute calibration of displacement amplitudes. 

In order to probe mechanical response, we take advantage of the optomechanical coupling 
and sideband resolved limit {ojm S> k) in an optomechanically induced transparency (OMIT) 
setup. In OMIT setup, two microwave fields are injected inside the cavity. A strong drive 
held pdrive at lower mechanical sideband frequency Ud = u)c — oJm and a weak probe held 
Pprobe measures the cavity response by sweeping the probe tone in the vicinity of Uc- When 
the detuning between drive and probe fields O = Up — cud matches Um, the mechanical 
resonator experiences coherent radiation pressure force. Coherent response of the mechanical 
resonator to the radiation pressure force up-scatter the drive held exactly at Up leading to an 
interference with the original probe held measuring the cavity response. This phenomenon 
is called optomechanically induced transparency (OMIT)^^^ as described schematically in 
Fig. 2a. It is worth pointing out that unlike heterodyne mixing schemes with low frequency 
RF drive, the radiation pressure force drive eliminates the need to apply a dc gate voltage. 
Furthermore, while the strength of the probe tone allows to control the driving force on the 
mechanical resonator, independently the drive tone can be used to tune the dissipation in 
the mechanical resonator using the optomechanically backaction. 

For an overcoupled single port cavity, the interference between the probe held and up- 
converted held (drive being at lower motional sideband) leads to an absorption feature at 
low driving powers as shown schematically in Fig. 2(b). In the linear response limit, the 
resulting rehection coefficient of the cavity can be written as, S'ii(a;) = l — where 

= -i{n-Jm)+'yn ./2 is the susceptibility of the mechanical resonator, = _i^n-L)+K /2 
is the suspectibility of the cavity, y™, is the mechanical dissipation rate, g = go-y/rid is the 
many-photon optomechanical coupling strength, and Ud is the number of the drive photons. 
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In the limit k S> 25f S> 7m, the measurement of the optomechanical induced absorption 
(OMIA) allows to directly probe the responsivity of the mechanical resonator giving its 
amplitude and dissipation rate, thus making it a sensitive technique. 

For a red-sideband drive Ud = ujc — oJm, the minimum value of the reflection coefficient 
is given by much simplified expression | ~ 11 1 where optomechanical cooperativity C is 

defined as C = In the limit of no optomechanical coupling ((7 = 0), we recover 12^ —1| 

expression for minimum for a single port reflection cavity, which sets the base line of OMIA 
feature. The linewidth of absorption feature is given by (1 -|- C')7m, where the additional 
term Cjm originates from the backaction effects of drive photons and can be tuned by n^. 
Furthermore, amplitude of the mechanical resonator can also be cast into a convenient form, 
X = Xzpf (y^) a/^- It instructive to see that for low cooperativity (O < 1), the 

mechanical amplitude can be tuned by both the probe and drive tone as a; oc ^UdUp. On the 
other hand, in the limit O > 1, the mechanical amplitude is proportional to suggesting 
that an increase in drive field leads to optomechanical damping and hence a reduction in 
the mechanical amplitude. An increase in the probe field, however, in both cases drives the 
mechanical resonator harder and yields larger amplitude. 

In Fig 3, we probe the OMIA response in detail by varying the number of intracavity 
probe photons hence the driving force, while keeping the number of drive photons fixed 
at Ud = 2.5 X 10^ and 1.0 x 10®. At low number of probe photons, the OMIA feature is 
determined by the linear response of the mechanical resonator. As Up is increased further, 
the nonlinearity in the OMIA response becomes evident with a stiffening of the mechanical 
resonator (positive shift in the resonance frequency) and the shark-fin like Duffing response 
accompanied by hysteresis with respect to frequency sweep-direction. 

In addition to the clear Duffing response, with the exception of the bottom two curves, it 
can also be seen that the OMIA dip on the non-linear regime becomes deeper. Qualitatively, 
the observation of a deeper OMIA dip when Up is increased can be understood from a 
reduction of the mechanical damping rate as the resonator is driven to larger amplitudes. 
Such a decreased mechanical damping rate would give a larger cooperativity and thus a 
deeper OMIA dip. In the last two curves, the cooperativity is continuing to increase, but the 
OMIA dip becomes less deep as the cavity has now crossed over to an effective undercoupled 
regime (see supplementary info of rePfor more details). In addition to the deeper OMIA dip 
that is suggestive of a decreased mechanical damping at higher drive forces. Fig 3 also shows 
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additional features. Comparing panels (a) and (b), smaller and larger Ud respectively, the 
mechanical linewidth in panel (b) is signihcantly larger. This is a consequence of increased 
optomechanical damping, which also explains the absence of hysteresis and shows only a 
transition to a Duffing response at higher power^. Finally, in panel (b), at the highest 
drive forces, we also observe an instability in the response in the form of a spike in reverse 
frequency sweep. 

To gain quantitative insight into these observations, we perform numerical hts on the data 
shown in Fig. 3. The nonlinear response can be primarily captured by including a Duffing 
term ax^ in the restoring force of the mechanical resonator^. Following ReP^, we perform 
numerical hts to extract the linear mechanical dissipation rate ( 7 ™), mechanical amplitude 
(xq), and the Duffing parameter (a) for any given probe and drive power. The gray curve in 
Fig 4(a) is the numerical htted curve overlaid on top of the experimentally measured data 
allowing us to extract the Duffing parameter to be a = 2.3 x 10^® kgm“^s“^. Using the 
analytical expression for the onset of Duffing bifurcation point Uup = 0Jm+ f m 
get a = 2.5 x 10^^ kgm“^s“^, which is close to the result we get by performing numerical 
hts. Fig. 4(b, c) plot the linear damping rate with the amplitude of the resonator extracted 
by performing numerical hts on datasets shown in Fig. 3. At low amplitude, we observe 
mechanical damping rates 7 ^ of 2ti x 700 Hz {Qm = 51760) for Ud = 2.5 x 10^, while 
for higher amplitude, the damping rate drops to 27r x 410 Hz {Qm = 88373). At large 
number of probe photons, the nonlinear dynamics of the OMIA becomes far more complex. 
The model with Duffing term in the restoring force still captures the response except the 
instability (sharp absorption feature in Fig. 3(b)) in the reverse frequency sweep. 

The decrease in observed damping rate at higher amplitudes suggests the presence of a 
negative non-linear damping term term in the equation of motion of the mechanical 

resonatoi^. As this negative nonlinear damping occurs also at low cooperativities, and 
as it is not seen in the theoretical calculations treating the optomechanical nonlinear re¬ 
sponse, we do not believe that it is an optomechanical effect, but instead intrinsic to the 
graphene resonator. There has been also observations of nonlinear damping in nanomechan¬ 
ical resonator^ and carbon based resonator^. One possible source of negative non-linear 
damping is the saturation of two-level-systems coupled to the mechanical resonatoi^^SMl At 
low drive powers, these two-level systems can absorb energy from the mechanical resonator, 
increasing the mechanical damping rate. At higher powers, the two-level systems become 
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saturated, and the damping rate goes down. Such an process was suggested as an expla¬ 
nation of power-dependent attenuation losses in glassed^^®, and also was used to describe 
power-depended dielectric losses in superconducting electrical resonator^. For such a sat¬ 
uration result in nonlinear damping effects, the level spacing of the TLSs should be larger 
than the bath temperature. In order for TLSs to describe the negative nonlinear damping 
observed here, the coupling between the TLSs and the mechanical resonator would have to 
be non-resonant, mediated by strong higher order processes. 

In conclusion, we examined the nonlinear dynamics of a graphene resonator coupled to a 
superconducting microwave cavity. In linear response limit, optomechanically-induced trans¬ 
parency measurements easily allows us to extract linear damping rate and peak amplitude. 
At moderate driving force when response becomes nonlinear, we perform numerical hts by 
including a Duffing term in the mechanical restoring force and hnd a = 2.3 x 10^® kgm“^s“^. 
Increasing the driving force further, the OMIA response becomes complex and it is no longer 
captured by the Duffing term. At these large amplitudes, higher order nonlinearities start 
becoming relevant and make the mechanical damping rate to appear low at larger ampli¬ 
tudes, where we observe a qualitatively new phenomena of negative nonlinear damping in a 
mechanical resonator. 


ACKNOWLEDGMENTS 

The authors would like to thank Andres Castellanos-Gomez, Sal Bosman, and Ben Schnei¬ 
der for their help during device fabrication and low-temperature measurements. The work 
was supported by the Dutch science foundation NWO/FOM and EC-FET Graphene Flag¬ 
ship. 


* g.a.steele@tudelft.nl 

^ J. S. Bunch, A. M. v. d. Zande, S. S. Verbridge, I. W. Frank, D. M. Tanenbaum, J. M. Parpia, 
H. G. Craighead, and P. L. McEuen, Science 315, 490 (2007), ISSN 0036-8075, 1095-9203. 

^ C. Chen, S. Rosenblatt, K. I. Bolotin, W. Kalb, P. Kim, I. Kymissis, H. L. Stormer, T. F. Heinz, 
and J. Hone, Nat Nano 4, 861 (2009), ISSN 1748-3387. 


6 




^ V. Singh, S. Sengupta, H. S. Solanki, R. Dhall, A. Allain, S. Dhara, P. Pant, and M. M. 

Deshmukh, Nanotechnology 21, 165204 (2010), ISSN 0957-4484. 

^ C. Chen, S. Lee, V. V. Deshpande, G.-H. Lee, M. Lekas, K. Shepard, and J. Hone, Nature 
Nanotechnology 8, 923 (2013), ISSN 1748-3387. 

® V. Singh, S. J. Bosnian, B. H. Schneider, Y. M. Blanter, A. Castellanos-Gomez, and G. A. 

Steele, Nature Nanotechnology 9, 820 (2014), ISSN 1748-3387. 

® P. Weber, J. Gttinger, 1. Tsioutsios, D. E. Chang, and A. Bachtold, Nano Letters 14, 2854 
(2014), ISSN 1530-6984. 

^ X. Song, M. Oksanen, J. Li, P. Hakonen, and M. Sillanpaa, Physical Review Letters 113, 027404 
(2014). 

® R. M. Cole, G. A. Brawley, V. P. Adiga, R. De Alba, J. M. Parpia, B. Ilic, H. G. Graighead, 
and W. P. Bowen, Phys. Rev. Applied 3, 024004 (2015). 

® A. Eichler, J. Moser, J. Ghaste, M. Zdrojek, 1. Wilson-Rae, and A. Bachtold, Nature Nanotech¬ 
nology 6, 339 (2011), ISSN 1748-3387. 

X. Song, M. Oksanen, M. A. Sillanpaa, H. G. Craighead, J. M. Parpia, and P. J. Hakonen, 
Nano Letters 12, 198 (2012), ISSN 1530-6984. 

A. Voje, J. M. Kinaret, and A. Isacsson, Physical Review B 85, 205415 (2012). 

X. Zhou, P. Hocke, A. Schliesser, A. Marx, H. Huebl, R. Gross, and T. J. Kippenberg, Nature 
Physics 9, 179 (2013), ISSN 1745-2473. 

M. Imboden and P. Mohanty, Physics Reports 534, 89 (2014), ISSN 0370-1573. 

J. D. Teufel, T. Donner, M. A. Gastellanos-Beltran, J. W. Harlow, and K. W. Lehnert, Nature 
Nanotechnology 4, 820 (2009), ISSN 1748-3387. 

G. Anetsberger, E. Gavartin, O. Arcizet, Q. P. Unterreithmeier, E. M. Weig, M. L. Gorodetsky, 
J. P. Kotthaus, and T. J. Kippenberg, Physical Review A 82, 061804 (2010). 

D. J. Wilson, V. Sudhir, N. Piro, R. Schilling, A. Ghadimi, and T. J. Kippenberg, 
arXiv; 1410.6191 [quant-ph] (2014), arXiv; 1410.6191. 

V. Singh, B. H. Schneider, S. J. Bosnian, E. P. J. Merkx, and G. A. Steele, Applied Physics 
Letters 105, 222601 (2014), ISSN 0003-6951, 1077-3118. 

G. S. Agarwal and S. Huang, Physical Review A 81, 041803 (2010). 

S. Weis, R. Rivire, S. Delglise, E. Gavartin, O. Arcizet, A. Schliesser, and T. J. Kippenberg, 
Science 330, 1520 (2010), ISSN 0036-8075, 1095-9203. 


7 



R. Lifshitz and M. C. Cross, in Reviews of Nonlinear Dynamics and Complexity, edited by H. G. 
Schuster (Wiley-VCH Verlag GmbH & Co. KGaA, 2008), pp. 1-52, ISBN 9783527626359. 

O. Shevchuk, V. Singh, G. A. Steele, and Y. M. Blanter, arXiv; 1507.06851 [cond-mat, 
physics:physics, physics:quant-ph] (2015), arXiv: 1507.06851. 

S. Zaitsev, R. Almog, O. Shtempluck, and E. Buks, in 2005 International Conference on MEMS, 
NANO and Smart Systems, 2005. Proceedings (2005), pp. 387-391. 

G. Zolfagharkhani, A. Gaidarzhy, S.-B. Shim, R. L. Badzey, and P. Mohanty, Phys. Rev. B 72, 
224101 (2005). 

F. Hoehne, Y. A. Pashkin, O. Astafiev, L. Faoro, L. B. Ioffe, Y. Nakamura, and J. S. Tsai, 
Phys. Rev. B 81, 184112 (2010). 

V. Narayanamurthi and R. O. Pohl, Reviews of Modern Physics 42, 201 (1970). 

P. W. Anderson, B. I. Halperin, and C. M. Varma, Philosophical Magazine 25, 1 (1972), ISSN 
0031-8086. 

W. Arnold, S. Hunklinger, S. Stein, and K. Dransfeld, Journal of Non-Crystalline Solids 14, 
192 (1974), ISSN 0022-3093. 

J. Gao, J. Zmuidzinas, B. A. Mazin, H. G. LeDuc, and P. K. Day, Applied Physics Letters 90, 
102507 (2007), ISSN 0003-6951, 1077-3118. 




FIG. 1. (a) A scanning electron micrograph of a multilayer graphene (10 nm thick) drum-shape 
resonator coupled to a superconducting microwave cavity (not shown here). Graphene resonator is 
suspended 150 nm above the bottom gate electrode, (b) Schematic diagram of the device: graphene 
resonator couples external microwave radiation to the cavity by forming a coupling capacitor. 
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FIG. 2. (a) Schematic showing the idea of radiation pressure driving. Due to the optomechanical 
coupling, driving the cavity with a strong tone near uid = ujc — and a weak probe tone ujp near 
Wc exerts a radiation pressure force on the mechanical resonator at ujm- The strength of radiation 
pressure force is controlled by the product of probe tone and the drive tone amplitudes, (b) Sketch 
of the cavity reflection coefficient in presence of a strong sideband drive. The optomechanical 
interaction produces an absorption feature in the cavity response, (c) A zoomed-in view of the 
OMIA feature showing the mechanical response in the linear regime. 
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FIG. 3. Forward (red) and reverse (cyan) frequency sweep measurement of OMIA feature show¬ 
ing mechanical response at various probe and drive powers. The probe photons are swept from 
Up = 2.5 X 10^ to 3.14 X 10® in 1 dB steps (top to bottom). Number of drive photons Ud is 
fixed at 2.5 x 10^ for panel (a) and 1.0 x 10® photons for panel (b). The evolution of nonlinear 
response accompanied by the hysteresis can be clearly seen as probe power is increased (top to 
bottom). Panel (b) shows instability points as sharp dips appearing at large probe power. For 
clarity, measurements in (a) and (b) are plotted with offsets of -30 dB and -9 dB, respectively. 
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FIG. 4. (a) Measurement of 5ii showing strong nonlinear response (red curve) together with 

numerically fitted curve (gray) for Ud = = 2.5 x 10^. (b, c) Extracted linear dissipation rate 7m 
plotted against mechanical amplitude for = 2.5 x 10^ and for 1.0 x 10®, respectively. 
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